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Abstract 



Bound states of heavy qq quarks are reviewed within the context of QCD. 

First of all, we consider the calculations which can be performed ab initio, which 

includes it with principál quantum number n up to four, bb states with n — 1 

and (to a lesser extent) bb with n = 2 and čc for n — 1. Among the results, 

f— j ' we report a very precise O(ai) evaluation of 6, c quark masses from quarkonium 

spectrum with a potential to two loops, a calculation of the decay T — > e + e~ and 
Q\ ■ a prediction for the splitting T — r/t- We then consider how, with the help of 

^^ ' reasonable assumptions, one can extend QCD calculations to other states of heavy 

^h , quarks. Finally, a few words are said on the treatment of light quark bound states. 
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1. Introduction 

In the present lectures we are going to review some aspects of the analysis of heavy quarkonia, it, čc and 
especially bb states. Before the advent of asymptotic freedom in 1973, hadronic interactions were analyzed 
with (among other methods) the help of the quark model, which incorporated a somewhat inconsistent set 
of semiphenomenological rules. An important role was played by bound statě calculations in the so-called 
constituent quark model, developed in the early sixties by, among others, Morpurgo, Dalitz and collaborators, 
and Oliver, Pene, Reynal and Le Yaouanc. In this model u, d, s quarks were given phcnomcnological masses 
of 300 — 500 MeV, and were bound by potentials: the harmonie oscillator potential bcing a popular choice 
because of its simplicity. Quite surprisingly, a large number of propertics of hadrons could bc reproduced in 
this way. 

After the advent of asymptotic freedom, and with it of a consistent ficld theory of strong interactions, 
it was possible to reformulate the quark model in terms of QCD. Thus, De Rújula, Georgi and GlashowW 
showed that taking into account relativistic corrections and colour algebra one could calculate the spectrum 
of the then known hadrons, including in particular such properties as the N — A mass differcnec, and even 
the S° — A splitting, something that had defied previous, non-QCD analyses. They were also able to predict 
qualitative features of the charmonium spectrum. 

Nowadays we expect more from QCD, at least for heavy quarks. The reason is that there, and to 
leading order in (v 2 ) ((v 2 ) the average velocity of the quarks), the interaction is equivalcnt to a potential. 
At very short distances this potential has to be of the coulombic type, 

-^. (1.1) 

r 

Even in the static limit, in QCD, (1.1) will be modified by radiative corrections; but these should bc 
of the form of a function of r. In fact, also at long distances one expects, in the nonrelativistic (NR) limit, a 
local potential with the form of a function U(r), although it will not be of coulombic type. The reason for 
this is galilcan invariance, that will hold in the NR limit. By virtue of it, we must háve that the derivative 
of the position Q, the velocity, should be proportional to the momentům: 

[H, Q] = ^Q = r^P, 

n nm 

with H the hamiltonian. If we define the interaction by H lnt = H — P 2 /2m and evaluate the commutator 
above, it follows that [iři n t, Q] = and so, because of a well- known theorem of von Neumann, H lnt must 
be a function of Q. A function which due to rotational invariance, and at least negleeting spin, may only 
depend on r: 

iři„t = U(r). 

Needless to say, relativistic corrections will in generál not be represcntable by local, momentum- 
independent potentials, as is the čase even in QED. In QCD one encounters QED-like corrections and 
idiosineratic QCD ones, in particular those associated with the complicated structure of the vacuum. Of these 
the most important are the effects involving the gluon condensate (a s : G 2 :), hrst studied in this context 
by Leutwyler and Voloshin™; the quark condensate also gives contributions but, for heavy quarkonium, 
subleading ones. 

We will consider in these lectures bound states of heavy quarks (and at the end say a few words about 
light quarks 1 ) in dcereasing order of tractability by rigorous QCD. First of all, we will consider situations 
where one has the inequalitics 

a<iž; \B n \-> A: m > A 

where a is the equivalent of the Bohr rádius, and B n are the equivalcnt of the Balmcr energies; m is the 
quark mass, and A ~ 300 MeV is the QCD parameter. Under these circumstances, nonperturbativc (NP) 
and confinement effects are expected to be small, and so are the radiative and relativistic corrections. All 
of them may therefore be treated as perturbations of a nonrelativistic, unconfined and leading order (in 
the QCD coupling a s ) systm. This will constitute the bulk of these lectures. Then we will diverge from 



Light quark bound states, and gluonic bound states, are reviewed in the companion lectures by Yu. Simonov. 



-HEAVY QUARKONIUM- 



the situations where one can effect ab initio QCD calculations, relying more and more on reasonable (but 
unproven) assumptions, and on models. 



2. Heavy quarks at short distances: pure QCD analysis 

For very heavy qq bound states the equivalent of thc Bohr rádius, a = 2/(mCVa s ), is much smaller than the 
confinement rádius, R ~ A . So we expect that, for lowest n states, with n the principál quantum numbcr, 
confinement may be neglected, or at least treated as a first order perturbation. In this čase the quarkonium 
systém is very similar to a familiar one, viz., positronium; so that methods analogous to thosc dcvcloped for 
thc last may be applied also to the study of the former. The NR potential may be obtaincd from perturbativc 
QCD; notě that, unlike for positronium, and because of the zero mass of the gluons, radiative corrections 
are present even in the static limit. Thc lcading piece is given by the tree level nonrclativistic amplitudě and 
we then include highcr cffects as perturbations: we proceed in steps. 

We will also follow the method of equivalent potentials, advocated for QCD by Gupta and collaborators^J :| 
we find this method the more transparent one. Othcr, equivalent methods, based on thc Bethe-Salpeter 
equation or effective lagrangians may be found in the literatuře ". 

In the method of equivalent potentials one profits from the fact that, in the NR limit, the potential 
is given by the Fourier transform of the scattering amplitudě, in the Born approximation: 



T 



A R 



'(p 



p) = 



1 



d re 



ir(p-p') 



V(r) 



(2.1) 



where p, p' are the initial and hnal momenta in the center of mass reference systém. On the other hand, 
Tm ™ may be calculated as the nonrelativistic limit of the relativistic scattering amplitudě: 



NR 



Tiborn 
1 NR 



= lim 



1 



-.F{p 1 +p 2 ^p' l +p' 2 )- 



formally the NR limit is equivalent to taking the limit of infinite quark masses, keeping the thrce-momenta 
fixed. One can thus calculate the Born approximation to F, F BoTn using thc familiar Feynman rules (actually 
at tree level) for qq scattering, také the NR limit and hence obtain T Bc £ n . From it, by inverting (2.1) one 
finds V. What is more, we can calculate corrections to V by including corrections (in particular, relativistic) 

iBorr 

NR 



to tP ™ 




FIGUŘE lA. A ladder of gluon exchanges. 

A few words will now be said about the connection of the method with Feynman graph calculations; 
more details may be found in Ch. 6 of ref. 5, which we roughly follow in these notes. Working in the 
strict nonrelativistic limit to avoid inessential complications, it is easy to check that solving the Schrodingcr 
equation with thc coulombic potential of Eq. (1.1) is equivalent to summing an infinite ladder of exchange 
graphs (Fig. 1a): indeed, both methods yield the samé S matrix. Including radiative corrections then alters 
the potential, as will be shown below. The corresponding Schrodingcr equation will bc equivalent to dressing 
the kernel Fn associated with the calculated radiative corrections to N loops with ladders, corresponding 
to the coulombic wave functions, as shown graphically in Fig. 1b where the kernel would be given, e.g. to 
two loops, by graphs depicted in Figs. 2 3 below. Thus the method of potentials is equivalent to a particular 
arrangement of thc summation of perturbation theory. 
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FIGURE 1b. Kernel, dressed with infinite sums of ladders. 



2.1. Step 1: coulombic Schródinger equation 

According to what we háve said, we proceed in steps. In this first step we consider a nonrelativistic tree-level 
kernel (Fig. 2). Here the scattering amplitudě is such that it generates a coulombic potential. So we get the 
Schródinger equation, for the energies and wave functions, 

#(<%(0) = £(0)^(0) 



1 



if(°) =2m A + y(°)(r), V w {r) 



Cf0í 8 



notě that for quarkonium the reduced mass is m/2. The quantities m, a s however are as yet undefined; 
only when including radiative corrections will we be able to give them a precise meaning. This equation, 
formally identical to that of an hydrogen-like atom, can be solved exactly, and will be our starting point in 
the calculations. 




FIGUŘE 2. One- 
gluon exchange. 



2.2. Step 2: relativistic corrections 

The relativistic corrections are identical to those for positronium, known since ancient times. They are 
found by considering still the tree level scattering amplitudě corresponding to the diagram of Fig. 2, but 
kecping now the terms of order 1/c 2 , c the speed of light; the details of the derivation may be found in 
textbooks on relativistic quantum mechanicsl 6 !. Adding also the correction to the kinetic energies, 



V r 



A/2m - A 2 /8m 3 



onc finds the hamiltonian, 

^«tree»=- ff( ° ) +^ ) ( 2 - 2a ) 

where the superseript zero in V^°' indicates that the potential is still obtained from a tree level (zero 
loop) amplitudě. The relativistic corrections, which are to be treated as first order perturbations to the 
unperturbed equation (2.1), read 



v: 



(«) 



rcl 



v 



(0) 



V, 



(") 



v, 



(0) 



LS 



v; 



(0) 



hf 



(2.2b) 
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The various pieces, spin- independent (in which we also include the correction to the kinetic energy), tensor, 

LS and hyperfine, are 

V {Q) = 1 A 2 CpOLs 1 

si;rcl 4m 3 m 2 r ' 

(0) _ C F a s 1_ 

tens _ 4 m 2 j.3* 512 ' 

o!™ t 2.2c 

^0) ^3^ 1^ 



(0) _ 47rC F q, 2 
Vhí ~ 3m 2 ^ Ó(rj - 
Here L is the orbital angular momentům operátor, S the total spin operátor, and S12 the tensor operátor: 



L = -irxV, S = ^l±^, S 12 = 2Vte-fc]S i S, 



The wave functions are assumed to háve spinor components, and the Pauli matrices a a act on spinor x(A a ), 
a = 1, 2. 

2.3. Step 3: radiative corrections 

Bcforc discussing the radiative corrections, a matter has to be settled hrst, which is that of the meaning of 
the mass in the Schródingcr cquation. We háve defined the potential by assuming that it vanishes at infinity; 
otherwise, we háve the ambiguity of an arbitrary constant. Thus, we must interpret the mass as the mass at 
long distances, i. e., on the mass shell. Now, both these requirements are not rigorously valid since quarks 
are confined in a region of rádius R ~ A ; but we can work with them to the extent that R is much larger 
than the region where the movement of the quarks takés pláce, a — 2/mCFCt s - 

With this requirement we define m, called the mass shell, or pole mass, to be such that, in pertur- 
bation theory, 

S P .t.W =m)- 1 =0. (2.3) 

The relation of this mass with the MS mass, m was found by Coquereaux and Tarrach^ to one loop and by 
Gray et alJ 8 ' to two loops. After correcting a misprint of the last reference one finds 

2' 
m(m 2 ) { 1 + ^ f ^ sv "" ; +(K- 2C F ) 



(2.4a) 



where, denoting by n/ the number of quark flavours with mass less than or equal to 



ni. 



-1 






..., .mJ' (2.4b) 




FIGUŘE 3. Some radiative corrections. 

All other quantities, however, are renormalized in the MS scheme. So, a s (/x 2 ) will be the MS coupling 
at the scale /i 2 , etc. These radiative corrections (some of which are shown in Fig. 3) háve been evaluated by 
a number of people. Those to the spin-independent part of the potential, in the strict static approximation, 
were hrst calculated by Fischler and Billoirel 9 ! to one loop and by Peter and SchróderI 10 ! (who checked and 
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corrected a mistake of Peter's calculation) to two loops. Because of the zero mass of the gluons, corrections 
nonanalytic in the average velocity occur. These depend on \v\ and are of a size comparable to that of the 
two loop static corrections; notě that, in the coulombic approximation, (\v\) = Cfol s . They were calculated, 
together with 0(v 2 ) radiative corrections to one loop, in rcf. 11. Spin-dcpendcnt one loop corrections were 
evaluated in refs. 3, 12. 

Wc discuss in some detail the spin-independcnt part of the spectrum. To také into account all tcrms 
giving corrections of 0(a^) to the energy spectrum one writes the Hamiltonian as 



H = ií (0) + H-í 
where H^ ' may, and will, be solved exactly and contains all the coulombic pieces of the interaction: 

m r 

7eA)\ ctslp 2 ) 



(2.5a) 



a a (/x ) =a s (/z ) { 1+ ai + 



7e 



*a +£)+(£ 



TV 



2 \ A) 



02 



H\ is to be considered as a perturbation. Its form is, 



#1 = 14i.;rel + V± L) + V^ ] + V {LL ^> + V s . Icl + V h{ , 



(2.5b) 



(2.5c) 
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s.rcl 
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47T 2 
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47r(7Fa: í 

3to 2 



s(s + l)í(r). 



(2.5d) 



^s.i.;rei is the spin- independent piccc of V rc \ in Eq. (2.2c); V s , re i is a one loop velocity-dependent correction 
and all other ones are one and two loop static corrections with the exception of the last term, representing 
hyperfinc splitting. Although we are considering the spin-indcpcndent interaction, to the precision we are 
working one needs to differentiate between the masses of vector and pseudoscalar states, hence the presence 
of this piece. cti was calculated in rcf. 9, ct2 in ref. 10 and b\ and many of the rest of the terms in ref. 11; 
all these constants are given in the Appendix. Notě that, to the precision we are working, a s has to be 
evaluated to three loops 



a s (Q 2 



1 



L 



47T 

= logQ 2 M 2 



& log L , 01 log 2 L-pl log L + f3 2 p -Ř 



P 2 L 



fáL 2 



Some bookkceping is necessary to identify which terms to include for a given order of accuracy in 
the calculation, e.g. in the evaluation of the energy levels. The pieces given here will provide a calculation 
accurate to order a*. All terms in H\ are to be treated as hrst order perturbations of H^\ except for 
the term V^ , which has to be evaluated to second order. Thus it produccs, in addition to the first order 
contribution, 






/3 C 2 a 2 (/x 2 )5 s (/x 2 
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log 2 



1 -7E 



(2.6a) 
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the second-order energy shift, for the ground statě, I 13 l 



S^IdEw = -m 



PoCj 



F^s 



4tt 2 



ATolog^+^logf 



No 



(2.6b) 



the constants iV are given in the Appendix. 

The hrst order contributions of the othcr V's are easily evaluated using the formulas of ref. 11. One 
finds 

V 1 



E^ ' = 2ra — m- 



y F^s 

4n 2 



The label "p.t." in Ef tl ' indicates that we háve as yet only ušed results deduced from perturbation theory; 
the full expression would be 

E n l 

with StqpE n i given below, Eq. (2.10). 



£ P t. 

ni 



5NpE rl l, 



The 5y E n i are, with a as above, 
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(2.8a) 



1-2) 



2r{n-l) 



- y 

i) z - 



r(2/ + 2 + j) 



r(n + l + í) j^ Q j\(n-l-j 
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CVM 2 



m n 3 (2Z + l)a 2 



Here and from now on we háve defincd a = 2/(CFina s ). We recall that constants are collcctcd in the 
Appendix. For the masses of the vector states (T, Y', T"; J/ip, ip' , . . .) one has to add the hyperfine shift, 
at tree level, 



r(l) 



E n 



SsíS, 



w 



8Cfcc s 
3n 3 m 2 a 3 



(2.8b) 



The value of the contributions of Vf were given above, Eq. (2.6) 



(L) 



2.4. Step 4: nonperturbative corrections 

The leading nonperturbative (NP) corrections can be shown to be those associated with the contribution 
of the gluon condensate. Physically they may be understood as follows. We consider that the quarks move 
in a medium, the QCD vacuum, which is full of soft gluons (Fig. 4) that we represent by their ficld strength 
operators, G c „ v (x). When a <C R, we may assume that the confinement size is infinite and, moreover, that 
one can neglect the fluctuations of the G'l l/ (x) in the region of size a in which the quarks move. So we 
approximate the effect of the motion in the gluon soup by introducing an interaction, which in the static 
limit will be of dipóle type, of the quarks with a constant gluonic ficld: 

H NP = - 9n G c 0l (0)ť = -grS c ť. 

Because of Lorentz invariance of the vacuum we assume that {G c ) — 0, but (a s : G 2 :) ^ 0. For dimensional 
reasons, this will give the leading NP contribution to the spin-independent energy shifts. Applying thus 
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FIGURE 4. The region where the qq pair 
move inside the confinement region. 



straightforward second order perturbation theory we havc 2 



S^pEni 



<%,H £ 



1 



Jí(«) - K 



(0) 



Be&X 



There are a few points to clarify regarding this equation. First of all, we may, since we average over directions, 
neglect the magnetic quantum number M . Secondly we háve ušed in the denominator the octet Hamiltonian, 

ff(8) = A + — — — . 

m 2N C r 

This happens because the perturbed statě, Hg^^) = — gr£ a t a \ x P^) , is manifestly an octet one, as £ a 
creates a gluon on top of the singlet \W^). 

Next, we háve to relate expectation values of products £ ... £ to the gluon condcnsate. For this, 
hrst write the gluon radiation Hamiltonian as 

1 



ífrad = 



8tt 



d d r : £ z 



B 2 



with sums over omitted colour indices undcrstood. Its expectation value in the physical vacuum should 
vanish, so we conclude 



(vac 



d J r : £ 2 : I vac) = -(vac I / d J r : B 



vac) 



We assume the field intensities to be constant, so we may replace the integrals by the volume timcs the 
integrands at x = 0. Canceling thcn the volume and recalling that G 2 = —2(£ 2 — B 2 ), we find 

(vac| : £ 2 : |vac) = -±(vac| : G 2 (0) : |vac). 

Finally, using Lorentz and colour invariance of the physical vacuum, 



4wa s dij6 a b 



g 2 (:£l(0)£í(0):)= " 



(D-Í)(N 2 -1) 



££:)=-^í(a s G 2 ) 



With this we get 



SNpE n i = U nh ír^t a £ a (0)j - 1 (Q) H £ ír^í 6 £ fc (0) j W n 

= ^TF^E {^un --A4 



-M 0) 



(2.9) 



rSf n 



To finish the calculation we háve to invert the operátor (— m _1 A + a s /6r). For the simple čase above, 
the method may be found in ref. 2; in more complicated situations, cf. ref. 14. So we finally obtain the 
nonperturbative pieces of the energy shifts, which are of the formel, 



$np E ní = in 



Trn 6 e n i(a s ■ G 2 
(mCFCts) 4 



(2.10) 



Projectors over the subspace orthogonal to íř^ '), that we do not write explicitely, are understood. 
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whcrc thc numbcrs e n i are of order unity, eio ~ 1.5. The evaluations for the spin-dcpcndcnt shifts may be 
found in ref. 14 (with a minor correction in ref. 15) and the contributions of higher order operators has been 
considered in refs. 16, 17. Notě that, as already remarked by Leutwylerl 2 !, one cannot derive (2.10) from a 
local potential; but, for the lowest states, the effect may be approximated by a cubic one, 

Vbluoncond.M ~ A 4 r 3 . 

As promised, the correction (2.10) is relativistic in that it is of order I/to 4 ; but the coefhcient is very 
large because of the high powers of n, ctj 1 . Thc rcason for these powers of a s and n can be understood easily. 
Two come from the energy denominators, and four from the expectation value {riTj) n i ~ n 4 /{mCFa s ) 2 . 

The right hand sidc of (2.10) grows as the sixth power of the radiál quantum number, n. It is in fact 
this very fast growth with n that leads to the breakdown of the method as soon as n exceeds, or in somc 
cases even equals, the value 2 for čc, 66; only for it can one go to n ~ 5. 

The NP corrections to the wave function may be obtained with methods similar to those employed 
to evaluate S^pE. For n = 1, l = 0, we háve 

#io(r) - (1 + ÍNp(r)) *io(r), (2.11a) 

where the NP correction is 

X (*.\ /2968 104 2 52 „3 1 „41 71"(q s G ) 2r 

Wr) = IW ~ 42šP - T275P - 225P ) m 4(c F a a ) 6 ' P ^~' ^ ' 

It turns out that the coefficient of the correction is larger than for the energy shifts, both in powers of aj 1 
and of n: the effects of confmement are larger for the wave function than for the energy levels. 

The contribution of some higher dimensional operators has been estimated by Pinedal 16 !; the cor- 
rections due to the finite size of the hadron is discussed in ref. 17, and will be briefiy reviewed later. 

Radiative and nonperturbative corrections to higher excited states may likewise be evaluated; as 
can be calculated the decay rates into photons or leptons, e.g., T — ► e + e~, rjb — * 77. We will present some 
of these results below. 

3. Results 

Let us summarize the results. The calculation is fully justified, in the sense that higher order corrections 
(both perturbative and NP) are smallcr than lower order ones for 66 with n = 1. The samé is partially 
true for the energy levels of the samé states with n = 2 and, for čc, for n = 1. For the wave functions of 
66, n > 2 and all čc states, and for the energy levels with higher values of n than the ones reported above, 
the calculation is meaningless as nominally subleading corrections overwhelm nominally leading ones. 

Bcfore presenting the results a few words háve to be said about the choice of the rcnormalization 
point, /i. As our Eqs. (2.6), (2.8) show, a natural value for this parameter is 

2 víiCfols 

Mo = — = , 

na n 

for states with the principál quantum number n, and this will be our choice. For states with n = 1 the 
results of the calculation will turn out to depend little on the value of u, provided it is reasonably close to 
(lo- Higher states are another matter; we will discuss our choices when we consider them. 
As input parameters we také the recent determinations' 18 ', 

A(n f = 4, thrce loops) = 0.283 ± 0.035 GeV [ a s {M 2 z ) ~ 0.117 ± 0.024 ] , 

and for the gluon condensate, very poorly known, the value 

(a s G 2 ) = 0.06 ±0.02 GeV 4 . 

(Notě that the slight difference between the results reported below and those of previous determinationsl 11 ' 13 ] 
are mostly due to the variation of the prefered value of A from 0.20 to 0.28 GeV.) 
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3.1. n = 1 states 



For bb one gets a precise determination of rrib, and less so of TOb(m 2 ) (pole and MS masses), a reliable 
prcdiction for the hyperfine splitting, and reasonable agreement with the experimental value of T — > e + e~; 
these will be discussed later. For čc a reasonably accurate value is also obtained for m c . For bb, and with A, 
(a s G 2 ) as given before and varying /i 2 around /it 2 , by 25% to estimate the systematic errors of the calculation 
one finds, from the T mass, the quark masses! 13 !: 

=5.065 ± 0.043 (A) T 0.005 ((a s G 2 )) 



m b 



+0.037 ( var y V 2 by 25%) ± 0.006 (othcr th. uncert. 



(3.1) 



m b (m 2 b ) =4.455 ± 0.012 (A) T 0.005 ((a s G 2 )) 

+o.olí ( var y V 2 by 25%) ± 0.006 (other th. uncert.). 

Notě that nu, is correct to 0(aj,([i 2 )), and fříř>(m 2 ) to 0{a 2 {m 2 )). The piece denoted by the expression 
"other th. uncert." in (3.1) refers to the error coming from highcr dimcnsional operators and some highcr 
order perturbative terms; it can be found discussed in refs. 13, 16. It is comfortably smaller than the errors 
due to the uncertainty on A, {a s G 2 }. 

The values of /z 2 ,, a s (/Zg), a s (^) are, respectively, 

tt 2 , = 7.86 GcV 2 , a s (/4) = 0.257 ,a s (/4) = 0.415. 

We see that a s is small, thus justifying the use of perturbation theory. Moreover, a/2 ~ (2.8 GcV) -1 <C A -1 : 
the quarks move well away from the confinement region. Finally, the binding energy is also considcrably 
larger than A. Thus wc find our approximations justified a posteriori. 

For čc, and with A(n/ = 3, thrcc loops) = 0.338 ± 0.037 and the mass of the J/tp as an input now, 



m r =1.936 



+0.059 
-0.068 



{A) T 0.014 ((a s G 2 )) 



+o.'ío6 (varying /i 2 by 25%) ± 0.014 (th. uncert.) 

.-'1 _1 K«/i+0.086 f ,n -,-n nl Q//„, r<2 



m c {m 2 c ) =1.564^;^ (A) T 0.013 ((a s G 2 )) 



(3.2) 



+[J;?^ (varying /z 2 by 25%) ± 0.013 (th. uncert.), 



and /íq = 2.871 GeV now. The errors, and the values of a s , a s inercase correspondingly and it follows that 
the calculation is much less reliable than for the bb čase, as the errors in (3.2) show. 

The values of the b quark masses reported here, e.g., Eq. (3.1), are slightly larger than those one 
finds with the sum rule method (see for example, refs. 19). It is not clear to the author why this oceurs. 
I suspect that the sum rule evaluations contain systematic uncertainties which are not under control; and 
indeed, the determinations are not very compatible one with another. Anyway, the diserepancies are not 
terribly large. 

3.2. n = 2 states 

For the states with n = 2, the energy levels can be evaluated using the values found for to& and taking 
now /i = l/a. However, since (as stated) radiative and nonperturbative corrections are large, the results are 
very sensitive to the value of \x chosen. For this reason it is more profitable to fít /i. This is the proceduře 
followed in ref. 14, from where the following table for the mass splittings of the states shown is taken: 



States 


Theory 


Experiment 


2 3 P 2 - 2 3 Pi 


21 ± 7 


21 ± 1 MeV 


2 3 Pi - 2 3 P 


29 ± 9 


32 ± 2 MeV 


2 3 S'i - 2 3 P 


181 ± 60 


123 ± 1 MeV 


2 3 S'i - l 3 S*i 


428 ±105 


563 ± 0.4 MeV 


2 3 P - 2 1 P 1 


1.5 ± 1 


- 



Here we use standard spectroscopic notation; the common, fitted value of n is somewhat above 1 GeV, and 
the errors given are only those generated by the errors in A, (a s G 2 ) given above. The overall agreement of 
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thcory and experiment is noteworthy, particularly considering that the value of the single free parameter, /z, 
is the samé for all states. However, the situation is less satisfactory than what a cursory glance to the Table 
might suggest: both perturbative and nonperturbative corrections are large, and the results are somewhat 
unstable. 

The wave functions for states with n = 2 present such large errors that the calculation using the 
methods described up till now become meaningless for them. 

3.3. Spin-dependent shifts and leptonic decay rates 

The evaluation of spin-dependent shifts, and decay rates follow patterns similar to thosc of the spin- 
independent energy shift evaluations. The expressions one finds are' 11 ' 14 ], 



M{V)~M{-q) = m 



Po 

2 



log T -l 



^(logC F a s 



1) 



q s 

Ti 



[1 + <^wf + t> N P 

1 Tr(a s G 2 



116 
8 704 



(3.3a) 



r(v - e+e-) = r(°) x [1 + Syri + S NP } 2 (1 + 6 lad ), 

2 



r (0) = 2 
4Cfcx 



Qb^QED 



M(V) 



(mC F a s (/j, 2 )y 



(3.3b) 



J rad 



3/? (, dfj, \ a s 



^NP 



1 [ 270 459 

2 L 108 800 



1838 781 1 
2 890 000 J 



Tr(a s G 2 



Here V = T, J/tp. The corrections are fairly large, particularly the radiative correctionl 20 ! <5 ra d. Because of 
this the calculation is less reliable than what one would háve cxpcctcd for 66, and fails completely for čc. 
With the values of rrib found before, one has the numerical results, 

M{r) - M{rfi = 53.3 ± 5.3 (A) ± 5.3 ({a s G 2 }) ± 10 {p 2 = 7.859 ± 25%) 



(3.4) 



and 



r(T -► e+e-) = 1.143 ± 0.11 {A) ± 0.11 ((a s G 2 )) ± 0.24 (/j 2 = 7.859 ± 25%). (3.5) 

Higher order NP corrections due to some higher dimensional operators are also known for the decay 
rate (see ref. 16). They would produce a shift in the decay rate of ~ 0.11 keV, smaller than the contribution 
of {ce s G 2 } or the uncertainty caused by e.g. varying \i around /iq- We do not include either in the evaluation 
or the error estimate. 

The calculated value for the decay is in reasonable agreement with the cxpcrimental figuře, 



^ exp. \-L 



1.320 ± 0.04 keV. 



4. Heavy quarkonia at long distances. Connection 
between the long and short distance regimes 

Herc we consider bound states of heavy quarks at long distances. This certainly includes čc with n > 1 and 
66 with n > 2; n = 1 for the first and n — 2 for the second are somewhat marginal. As stated in the previous 
section, perturbative QCD supplcmented with leading NP effects fails now; but, fortunately, and since the 
average velocity of bound states deereases with inereasing n, we expect the dynamics to be governed by a 
potential: our task is to determine it. This has been considered by a number of pcoplcl 17,21 ' 22 !. Here we will 
follow the Dosch Simonov method, in the version of ref. 17, which will allow us to establish connection with 
the short distance analysis of the previous section. 
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The potential, that we denote by U(r), is expected to exhibit a number of features. First of all, 
it should behave as ar at long distances, as follows from e.g. the lattice calculations. Secondly, it should 
contain a coulombic piece, so we write 



U(r) = -- + U NP (r), (4-1) 

r 



and, at short distances, one idcntifies k — CpOis + radiative corrections. 

To find this potential we consider the Green's function in terms of the Wilson loop, working directly 
in the nonrclativistic approximation, and for large time T: for a qq pair: 



G(x,x;y,y) = / VzVž C -^ Ka+Ka \W{C)), 

T ( 4 - 2 ) 

(W(C)) = [ VBe$» dt(L ™ t+L ™ d) . 



Notě that we treat the quarks in the nonrclativistic quantum mcchanical formalism, appropriate because 
of thcir nonrclativistic motion. Thus, Ko,Kq & r< 3 time integrals of the kinctic cncrgics (nonrclativistic 
lagrangians) of quark and antiquark, 

i-T „, r T 



Ko = ^J dž(í) 2 , K Q = ^f dÉ(í) 2 . 



* JO 

2 



However, the gluons are treated fully field-theoretically. So the radiation lagrangian is L ra( i = — i J d 3 rG 
The Wilson loop operátor corresponds to the contour C enclosing the q, q paths from time to time T. It 
should include path-ordered parallel transporters for the initial and finál states, Š>(x,x), <P(y,y) with e.g. in 
matrix notation 



<f>(x,x) — Pexpig / 

J X 



áz^t a B a Jz) 



which we do not write explicitly. 

To také into account the nonpcrturbative character of the interaction it is convenient to work in the 
background gauge formalism and write B^ = b^ + a^ where the a^ represent the quantum fluctuations and 
b^ is a background ficld. This is constructed such that the vacuum expectation value of the Wick ordered 
products of the a^, and of the mixed a^, b^ products vanish. Therefore, we may express the gluon correlator 
in terms of 6 M only: 

(:G(x)G(y):)^(:G b (x)G b (y):), 
G btllv = d^b v - d u b^ + gb^ x b v . 
Expanding in powers of the background ficld b^ we may write the Wilson loop average as 

(W(C)) = í VaPe-fc dz " a » 

U/ ] ^ Va J dz M í áz' v P<P a (z, z')b^z)P<P a (z' , z)K{z') + ... (43) 

=W + W 2 + ... 

and the transportér <P a is constructed with only the quantum field a. For the first term, Wo, the cluster 
expansion gives 

Wq = Z exp ((f2 + highcr ordcrs) , 

C ^ f T étí T éť 1 + ZZ ' 



4tt 2 Jo 7 r 2 + (í-ť 

r T 

-1 



¥>2 

**'" J '■ 

^Cpasr- 1 I át + 0(v 2 
Jo 

i.e., the coulombic piece of the potential. (Z is a constant that, in particular, includcs rcgularization) . 

We keep now only the next piece, Wi and neglect the W n , n > 2. At short distances, this is 

justified because highcr W n involve higher powers of the background fields. Thus the ensuing corrections 

- 11 - 
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are suppressed, on dimensional grounds, by powers of í/m n . For long distances, and although argumcnts 
háve been advanced for the dominance of the W2, we really háve a model, the so-called "stochastic vacuum 
model" . 

For details of the evaluation of this first nontrivial piece, W2, we refer to the lectures by Simonov. 3 
It produces a correction to the Green's function, SG, which in the static approximation is 



5G = -£ í d 3 r í dV í r t df3 l r\d0 
xG^(r(T),r)GP(r,r')G { c S \r',r(0)). 



G c ' are the singlet, octet coulombic Green's functions. The reason for the appearance of G^ 8 ' are similar 
to those for the appearance of H^ in (2.9). 

We may then také matrix elements between coulombic states, \nl), and identify the cnsuing cnergy 
shifts from the relation 

G = G { c ] +5G ~ G ( c\l-T6E nl ). 

T — >-oo 

We then find the basic equationt 17 ! , 

Mm = ± /"^ř ídf3d(3'Á( P )J2(nl\ne^-^\k(8)) 

J ( ] J (4.4) 

x-^ (fc(8)|ríe i "'(' í - 1 / 2 )' , |^). 

El' - E n -p 

(8) 

The states |fc(8)) are eigenstates of the octet Hamiltonian, with energy E) ; the E n are the coulombic 
energies. 4 Finally, A(p) is dcfincd in terms of the correlators, being the Fourier transform of 

A(x) = D(x) + £>i(x) + x 2 d 2 D 1 (x)/dx 2 

and 



( 3 2 :G 0í (x)G ,(0):) = ^ 



SijD(x) + XíXj 



d 2 D 1 

dx 2 



We may write, using Lorentz invariance, A(x) — A(x 2 /T 2 ), with T g the so-called correlation time. This will 
play an important role in what follows. 

We háve now two regimes. If \xt = Tg 1 3> \E n \ the velocity tends to zero, and the nonlocality also 
tcnds to zero as compared with the quark rotation period (which in the coulombic approximation would be 
1/l-EVtl). We can now neglect, in Eq. (4.4), both E n , E k ' as compared to po so, after some elaboration, we 
obtain the energy shifts 

5E n i ~ (nl\U^p\nl) 

with 

ir ( r f°° 
f/ Np(r)= — { / dXj dvD(\v) 



+ i j d\ 2 í dv [-2D(A, v) + D 1 (X, v)} }, 



(4.5) 



/o ./o 

-n/J „2 



D(X,u) =D(xq,k' s ), etc. 

At large r, and as this equation shows, we find U-^p(r) ~ ar. Here a can be related to T g and the gluon 
condensate if we assume a model for A. So, if e.g. we také an exponential ansatz for A{x), as in ref. 22, we 

find 

■k (a s : G 2 :) 
\á,t = — 7= i — 0-32 GeV . 



3V2 



a-i 



3 For technical reason, should be performed in euclidean space returning to Minkowski space at the end of the 
calculation. We omit these subtleties here. 

4 Here we subtract the rest energy, 2m, from both E^ 8 ' , E n . 
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For small r the limit of the expression in Eq. (4.5) leads to' 21 !, 

í7 N p(r) ~ c + cir 2 . 



(4.6) 



This is different from the behaviour expected from the Leutwyler-Voloshin analysis which gave a behaviour 
~ r 3 ; but one should understand that the present derivation holds for r — ► but still T" 1 ^> \E n \, i.e., in a 
situation other than that where the Leutwyler-Voloshin analysis is valid. 

It may be stated that the analysis based upon the potential (4.5) gives a very good description of 
heavy quarkonia states' 23 '. Notě, however, that the description is not the only one available on the market; 
othcrs, based e.g. on relativistic corrections to an assumcd lincar potential are given in the papers in ref. 24. 

We next get the matching between the two regimesl 17 ! . For this we turn to the opposite situation, 
viz., TV 1 <C \E n \. Now we may approximate A(x) ~ constant, so that A(j>) ~ 5i(p) and Eq. (4.4) becomes 



^NP-EnZ — 



:G 2 



18 



-(nl\i 



1 



#( 8 ) - E„ 



Ht 



il), 



(4.7) 



which coincidcs exactly with the results of the Leutwyler-Voloshin analysis^ 2 ! in the limit T g — > oo (/zt — > 0): 
cf. Eq. (2.9). In fact, Eq. (4.7) allows us to estimate the finite size corrections to the Leutwyler-Voloshin 
NP effects, which improves still the agreement between theory and experiment ^ 17 '. 

5. Further discussion of nonperturbative effects: 
Renormalons, and saturation 

In the previous sections we háve shown how QCD can give a very satisfactory account of the heavy quarko- 
nia spectra, particularly of the lowest lying states; an understanding based on perturbative calculations 
supplemented by NP ones, in particular those associated with the gluon condcnsate. Here we address some 
questions related to that. 

First, one may inquire about the connection of renormalons with nonperturbative effects. We return 
to the one-gluon exchange diagram, Fig. 2. If we dress the gluon propagátor with loops as in Fig. 5 then the 
corresponding potential, in momentům space, is 

V(k) = — — pQlog{k 2 /A 2y (5- 1 ) 

and we háve substituted the one-loop expression for a s (k 2 ). The expression (5.1) is undefmed for soft 
gluons, with k 2 ~ A 2 . As follows from the generál theory of singulár functions, the ambiguity is of the 
form cS(k 2 — A 2 ): upon Fouricr transformation this produces an ambiguity in the x-space potential of 
SV(r) = c[smAr]/r. At short distances we may expand this in powers of r and findl 25 ' 26 ! 



SV(r) ~C + C\r 2 + 



(5.2) 



The samé result may be obtained with the more traditional method of Borel transforms. The behaviour in 
Eq. (5.2) coincides with the short distance behaviour of the nonperturbative potential U^p(r) as determined 
in ref. 21, and Eq. (4.6) here. 




FIGUŘE 5. One- 
gluon exchange, 
dressed with loops. 
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FIGURE 6. Emission and absorption of a soft gluon collectively by a qq pair. 



The situation just described applies for states qq at short distances; but not so short that zero 
frcquency gluons cannot separate the bound statě. If this last is the čase, soft gluons do not resolve the 
qq pair and only see a dipóle. The basic diagram is no more that of Figs. 2, 5, but that of Fig. 6. The 
generated renormalon may then be seem 26 ! to correspond to the contribution of the gluon condensate in the 
Leutwyler-Voloshin mechanism. 

The matter of renormalons does not end here. If we calculatc the renormalization of the mass to 
one loop, and dress the propagátor with bubblcs, one also finds a renormalon contribution to an ambiguity 
in the mass oť 27 ' 

^renormalons ^ /^; 

and there is another nonperturbative correction to the hamiltonian related to the arbitrarincss in the origin 
of the energies. Indeed, we fix this origin by requiring the potential to vanish at infinity but, since the quarks 
are confined in a region of rádius R ~ A~ , "infinity" is equivalent to R, hence we get an ambiguity of order 
A ~ 1/R. 

The situation, however, is less confused than what one might think. In fact, and at least in the čase 
in which m\v\ is large compared to A, it can be shown that the linear and quadratic renormalons in the 
potential and pole mass cancel, leaving a yl 4 r 3 renormalon. This was to be expected on generál grounds (see 
e.g., ref. 26); a formal proof may be found in ref. 28. 

To make matters worse (or to improve them, depending on the viewpoint) we will also consider the 
possibility of saturation. We notě that the ambiguities we háve found are associated with small momenta 
or, equivalcntly, long distances. However, at least the singularities are clearly spurious. Indeed, not only 
the theory should be wcll defined but, becausc of confincment, long distances are nevěr attained: the theory 
possesses an internal infrared eut-off of the order of the confinement rádius, R ~ A . To try and implement 
it we consider again the gluon propagátor. To one loop it gets a correction involving the vacuum polarization 
tensor. Neglecting quarks this is, in x-space, given by an expression likc 



n${x,o) 



^9 Jabcja' de 



d 4 yi d 4 y 2 TBZ( yi )d^B ca ( yi )B%(y 2 )d„B el3 (y 2 )\0} 



We can také into account the long distance interactions by introducing a string between the field products 
at finitc distances. In matrix notation for the gluonic fields, B 11 — t a B£, this is implcmcnted by replacing 



B a { yi )B' 3 {y 2 ) -> B a ( yi )P (exp ij m WB^z)) B?(y 2 



The process may be described as "filling the loop" (see Fig. 7) by introducing all exchanges between 
the gluonic lineš there. If we furthermore replace the perturbative vacuum |0) by the nonperturbative one 
|vac), then a calculation similar to that made for the long distance potential for heavy quarks in Séct. 4 
yields a dressed propagátor 



DZ sscd (k) = D^(k) 



4li 



/3 log(M 2 + /c 2 )/yl 2! 
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FIGUŘE 7. 
gluon loop. 



"Filled in" 



and M 2 is related to the gluon condensate at finite distances, (G(x)G(0)) vac . 

This indicates a saturation property of the coupling constant at small momenta (long distances); 
the calculation in fact suggests that, at small momenta, the expression for the running coupling constant 
should be modified according to 

4-7T 4-7T 

*Ák 2 ) = -TT-TTT7-, - ^(fc 2 ) = 5-T-7I5 , ^ t Á2 - (5.3) 



P logk 2 /A 2 



/3 log(fc 2 + M 2 )/yl 2 ' 



It is certain that an expression such as (5.3) incorporates, to some extent, long distance properties of the 
QCD interaction. For example, if we také (5.3) with M = A in the tree level potential for heavy quarks, this 
becomes the Richardson potentiaP 9 ' 



yW(k) 



47rC F a s (k 2 ) 
k 2 



y(0),sat( k ) _ 



16n 2 C F 



/3 k 2 log(fc 2 + yl 2 )/yl 2 



When one has k 2 3> Á 2 , the short distance coulombic potential is, of course, recovered. For k 2 <C A 2 , 
however, 

167T 2 C F /1 2 



y(0),sat( k ) 



whose Fourier transform gives 



y(o),sat/ r \ 



k 2 «yl 2 /3 k 4 

(constant) x r, 



r>yl 



i.e., a linear potential. Indeed, a reasonably aceurate deseription of spin-independent splittings in quarkonia 
states is obtained with such a potential. Likewise, use of (5.3) with M — A provides a surprisingly good 
deseription of small- x deep inelastic scattering down to Q 2 ~ 0, as discussed in ref. 29; and these two cases 
are not unique. 

In špite of these successes, it should nevertheless be obvious that (5.3) can only bc of limited 
applicability. For example, consider the correlator of two currents in the spacelike region, II(Q 2 ). We 
know that in some cases such as the correlators of vector or axial currents for massless quarks, or that of 
pseudoscalar ones, one has 

n(Q 2 ) ~ 77pcrturbativc{l + 0((a s G 2 ))Q- 4 }, 
Q 2 ^oo 

whereas (5.3) would give a correction of order M 2 Q~ 2 . The Richardson potential is also a good example of 
the limitations of the uses of saturation, in particular in connection with the extent to which saturation really 
does (or does not) represent a reál, physical improvement, or merely the addition of a somewhat arbitrary 
new parameter. Indeed, the linear potential induced by saturation in the Richardson model is the fourth 
component of a Lorentz vector, while we know that the Wilson linear potential, as obtained, e. g., in the 
stochastic vacuum model or in lattice calculations, should be a Lorentz four-scalar: it thus follows that the 
linear potential obtained from saturation can be only of phenomenological use in some specific situations. 

It is not easy to draw a clear morale from all of this. One can try to eschew the problém by expressing 
observables in terms of observables (for example, r(T — ► e + e~) in terms of M(T)), hoping that this will 
reduce renormalon ambiguities, as some calculations seem to indicatel 27 ' 30 ! . This is the viewpoint adopted in 
the papers in ref. 31. Another possible attitude is the following. It is very likéry that the perturbative series 
in QCD are divergent; hence, different methods of summation lead to different results. This appears to be the 
čase for renormalons or saturation resummations. It is the author's belief that only if the summation method 
is rooted on solid physics it is likely to represent an improvement; otherwise, estimates of nonperturbative 
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cffccts become pure guesswork. In this respect, the method of taking into account thc nonpcrturbativc 
nature of the physical vacuum by considcring the effect of nonzero values for the correlators stands somc 
chance of being meaningful, as indccd phcnomenological calculations seem to indicate. 

6. Models 

6.1. The Constituent Quark Model 

We hrst discuss the constituent quark model. Here, we assume that the fact that quarks inside hadrons 
move through a medium made up of gluons and quark-antiquark pairs can, under certain circumstances, 
be represented by ascribing an cffcctivc mass, called thc constituent mass, even to light quarks. A possiblc 
way to connect this with a QCD analysis might be the following. 5 Consider the quark propagátor, in the 
physical, nonpcrturbative vacuum that we denote by |vac), 



Sij(x) = (vac\Tq i (x)q 1 (0)\ 



vac 



it is a gauge dependent object. We may derme an invariant propagátor by inserting a line integrál. In matrix 
notation, but working in Minkowski space for now, we thus write an effective propagátor as 

Sff(x) = -^(vac\Tq(0)Pcxp- [9 I" dy " B ^ y) g(aO|vac>. (6.1) 

N c 

We can interpret Sfj* as the propagátor describing a quark as it moves in the gluonic soup inside a hadron. 
In p-space, 



COff 



(p) = Já 4 xe i ^S^(x). 



Stf^)^jJ^-^---^-(q q )}; ((i.2) 



Evaluating the short distance limit with the OPE, we get the familiar lowest order expression 

? tf (x) x-o 6ij { 1^ x^o " M 

for simplicity we háve taken the quark to be massless. At long distances we evaluate Sfj(x) as follows. First, 
we go to Euclidean space. Thcn, and because we expect confincmcnt (and thus that the interaction grows 
at long distances), we calculate for large coupling, g — > oo. Finally, the quenched approximation is ušed. 

Under these circumstances, the evaluation of Sff (x) is identical to that of the Wilson loop. Undcr- 
lining Euclidean quantities, we then find 

Sf?(x) ~ ^e- 1/2 l2l ; 

J x — *oo 

a is the string tcnsion. In Minkowski space this becomcs 



Sg(x) ~ Šije-"' <-*\ (6.3) 



an expression which is very appcaling. According to it, the probability of a quark in the vacuum (inside a 
hadron) to propagate at a spacclike distance r = \J — x 2 decreases exponentially when r 3> cr -1 ' 2 ; but the 
quark may move frčely along a timelike or lightlike trajectory, where \/—x 2 is pure imaginary or zero. 
A simplc ansatz incorporating short and long distance behaviour is 

The corresponding p-space expression is then easily evaluated to be 



''■> l ' \J V" (a - p 2 - ioy/ 2 ) ~ n c (k -p 2 - ío) 5 / 2 



^)=^ i- , ; a in , 1/a - „z - a ™l„ ■ (6.5b) 



Other mechanisms for the generation/interpretation of a constituent mass may be found in the lectures by Yu. 
Simonov, and in ref. 32. 
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The expression (6.5) for the propagátor fulřills thc Bricmont Frohlicm 33 ! critcrion for confincmcnt 
and indeed exhibits many of the characteristics of the propagátor for a particle with nonzero effective mass. 
Thus, Slf(p) presents a cut starting at p 2 — K and, what is more interesting, it behaves for p — * like the 
propagátor for a massive particle: 

^)- 4 + ^¥---> (6-6) 

where the effective mass po is 

N C <J 2 
-ÓTT z (qq) 

It is curious that in the last expression the quark condensate appears in the denominator. 

The numerology works reasonably well. With the value po ~ 320 MeV obtaincd from phcnomcnolog- 
ical quark models, we can predict a 1 / 2 ~ 470 MeV, in reasonable agreement with lattice calculation results 
that give a 1 / 2 ~ 420 MeV. These nice features should, however, not makc onc forget the shortcomings of 
our calculation herc; (6.5) is to be considered no more than a phenomcnological expression. In fact, not 
only is the interpolation ušed somewhat arbitrary, but, because the expression for the propagátor only takés 
account of a certain class of gluon couplings, use of (6.5) into Fcynman diagrams may lead to violations of 
gauge invariance. Because of this it is probably better not to ask too much of the model and také, simply, 
the consequence of the existence of a universal mass that represents the inertia acquired by quarks due to 
their having to drag in their motion the gluon-quark soup, and which adds to mechanical quark masses. 
Concentrating on light quarks, we then assume masses 

m u (const) = m u + Ho, md(const) = rrid + po, m s (const) = m s + po- (6-7) 

The presence of the mass po breaks chiral invariance, and therefore pions and kaons (in particular) 
will be very poorly deseribed by the constituent quark model: for these particles we háve to use different 
methods. But one can use the constituent quark model to deseribe with success other hadrons (p, K* , S, 
A, nucleons, A, . . . ). 

To implcment the interactions among quarks, wc introduce two phenomcnological potentials: a 
confining potential, lincar in r, 

U CO ní(r) = Ar, A - a, (6.8a) 

and a coulombic-type interaction, 

tfcaui(r) = — , (6.8b) 

r 

together with corresponding QCD-type hyperfinc interactions. For quarks with indices i, j, we také 

U hyp (r) = -kJ2 —^— E *?*i ff <^. ( 6 - 8c ) 

and ti, o i act on the wave function of quark i. k may be connected with thc running coupling at, say, the 
reference momentům of 1 GeV: 

K~C F a s {l GeV 2 ). 

Because the model is in any čase not terribly precise, one at times replaces the linear potential by a quadratic 
potential, which can be solved explicitly. 
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6.2. The Bag Model 

In the bag model one neglects confincment, and replaces it by the boundary condition that quarks cannot 
leave a "bag" with size R ~ A^ 1 . There are a number of variants of the bag: the quantum mcchanical 
bag, the field theoretic bag, little bags, etcJ 34 ' 35 ' The simplest, of course, is the quantum mechanical bagl 34 '. 
There, as a first approximation, one considers quarks as free, subject only to the boundary condition that 
the wave function vanish for r = R. This bag is solved by writing the Dirac equation for free particles and 
imposing thcn the boundary condition. The simple model produces qualitativcly reasonable results for light 
quark states, with the exception of pionic and kaonic states: as was to be expected, because the presence 
of the bag breaks chiral invariance. It is also not clear to what extent the presence of the bag is a good 
simulation of the confincment mechanisms or, more gcncrally, nonperturbative cffccts: in fact it is not for 
heavy quarkonia. To see this, consider a systém of heavy quarks, qq, inside a spherical bag of rádius R. We 
také the interaction to be the coulombic one, —Cpa s /r, and impose the bag boundary condition. 

Let us denote by m, r to the reduced quark mass, m r = m/2 for quarkonium. For a statě with energy 
E, derme Ue = ^/Ry/(— E), and the variablc p = 2r/riE<i, where a = í/m r Cpa s and Ry = ^m r (CECt s ) 2 . 
For states with l — 0, the diffcrcntial equation obeyed by the radiál wave function We is 

P"(P)+ (- l) P'{P) + ^-^(p) = o, 
\P ) P 

V E {p) =(Const.) x c-p' 2 P{p). 

Moreover we háve the bag boundary condition P{L) = 0, where L = 2R/riE(i- 

We are interested in the solution of this for L>a. To find it, we proceed as follows. The regular 
solution is proportional to Kummcr's function, P{p) = M{\ — nE,2;p); the boundary condition then fixes 
tie- To sec this, consider the ground statě. Since in the limit L^oowc should recover the ordinary solution, 
with ue = 1, we writc tie = 1 + e and work to lowest order in e. Expanding, 

M(í -n E ,2: p) =M(-e, 2; p) ~ M(0, 2; p) + e^M{a, 2; p)\ a=Q 

1 v—v r(n) „ , . 

r (- € )Žť 1 nlr ( n + 2 ) 



v(p) 



^c 



^— í n(n + 1)! 

n— 1 v ' 

For this to háve a zero at p = L we must háve 

L n 



e- 1 = n(L) = J^ 



' n{n + 1)! 
At large L, rj(L) ~ e L /L 2 , hence e ~ L 2 cT L = (2R/a) 2 c~ 2R / a . For the energy, therefore, 

L — >oo 

E = -Ry + 2mlR 2 {C F a s ) 2 c- 2R ' a . 

In the čase of quarkonium, the nonperturbative shift induced in the ground statě by the presence of the bag 
of rádius R is thus 

S NP E = lm 3 R 2 (C F a s ) 2 e- RmCFa % 

totally different from what is found both at short distances (as caused by the gluon condensate, Eq. (2.10)) 
or at long distances, as given by a linear potential ar that yields something proportional to a/ma s . 
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Appendix: Constants 

We collect here, for ease of reference, some of the constants that appear in the radiative corrections to heavy 
quarkonium states. 

0o = 11- |n/;/3i = 102- f n f 

o _ 2857 5033 ^ , 325^2 



3ÍC A - 20T F n f C F - 2C A 
cii = ^7 — 1-47; oi = ~ —2.33; 



36 



" [W + ¥ C(3)] C A T F n s [f - 16C(3)] C F T F n f + f^}; 



4 L) = oiA, + ijSi + Í7E0O 



2 



5 = |(1 - log2)T F - lT F n f + il^_^£. 
^ = V>(l + f + n)-l . ^(„.o = 1^ (1 + z + n) ^ (1 +l + n) _ 2 - 

(n + l)\ ^ sl(s + l + l- n) 3 

x ' s—0 v J 




(n-l-iy. ^ (s + 2l + iy.(s + l + l~n) 3 

We remark that the definitions xof a^ and &i háve been swapped, in contrast to refs. 11, 13, 14, but 
in agreement with refs. 10. 
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